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The probabilities P±(to, i) that a growing Kardar-Parisi-Zhang interface remains above or below the 
mean height in the time interval (to,t) axe shown numerically to decay as P± ~ (to/t) ± with 6+ = 
1.18±0.08 and 0- = 1.64±0.08. Bounds on 9± are derived from the height autocorrelation function 
under the assumption of Gaussian statistics. The autocorrelation exponent A for a d-dimensional 
interface with roughness and dynamic exponents /3 and z is conjectured to be A = (3 + d/z. For a 
recently proposed discretization of the KPZ equation we find oscillatory persistence probabilities, 
indicating hidden temporal correlations. 



PACS numbers: 05.40. +j, 81.10.Aj, 02.50.-r 

The notion of persistence plays a central role in the 
assessment of random processes in science and everyday 
life. Intuitively, a process is persistent if it tends to main- 
tain its current trend. For example, if a continuous ran- 
dom process with zero mean is found at a positive value 
at time to, one may ask for the persistence probability 
P+(to,t) that it remains above zero up to time t > to] 
the probability P_ (to, t) for remaining negative is defined 
analogously. 

Here we consider the persistence of growing interfaces. 
The random variable of interest is the height h(x, t) above 
some fixed but arbitrary substrate site x, relative to the 
mean height h(t). At time t = the interface is flat, 
h(x,t) = 0, and we wish to characterize P±(to,i) for 
< to < t; in general positive and negative persis- 
tence probabilities differ, because the growth breaks the 
h — > —h symmetry. The dynamic scale invariance of in- 
terface growth jjj implies a scaling form for the temporal 
autocorrelation function, 

A(t, t') = ({h(x, t) - h(i))(h(x, t') - h(t'))) = t 2f3 A(t'/t), 

(1) 

where (3 > is the temporal roughness exponent; in par- 
ticular, the surface width W(t) — \J A(t, t) grows as t 13 . 
By the same token, the persistence probabilities P±(to, t) 
can depend only on the ratio of the two time arguments, 
and for t 3> to one expects a power law decay 



P±(t ,t) ~ (t /t) 



(2) 



defining a pair of persistence exponents 9±. 

A natural question concerns the relationship, if any, 
between the roughness exponent (3 and the persistence 
exponents 9±. If (3 is large, the height quickly deviates 
from the mean and is unlikely to return; thus, 9± should 
be small. Indeed, when the initial condition h(x, 0) is 
an interface with fully developed (stationary) roughness, 



the translational invariance in time can be exploited to 
derive the simple relationship [^|j3| 



= 1-/3, 



(3) 



where the superscript refers to persistence in the steady 
state. For the flat initial condition a detailed study Q 
of linear growth equations (for which 9 + = 0- = 6o) 
shows that #o still increases with decreasing /3, however 
the precise dependence is nontrivial and only accessible 
to perturbation theory and rigorous bounds. In fact the 
two exponents 9^ and (?o describe different asymptotic 
regimes of a single generalized persistence probability || . 

In the present paper we extend the work of Ref. Q to 
the generic nonlinear interface growth equation |l| 



dh 
~dt 



A 



vV 2 h+ -{Vh) 2 +r)(x,t) 



(4) 



introduced by Kardar, Parisi and Zhang (KPZ) Q. This 
entails several complications. First, we expect (and ver- 
ify numerically) different values for 9+ and Second, 
due to the nonlinearity of ([|), the height fluctuations are 
non-Gaussian and therefore the autocorrelation function 
(|l|), some properties of which are known [g— g| , does not 
even in principle fully determine the persistence proba- 
bility. The techniques developed in Q are used to derive 
bounds on the persistence exponent of a Gaussian process 
with the KPZ autocorrelation function. Our numerical 
estimates violate the bounds, indicating that the non- 
Gaussian nature of the interface fluctuations is impor- 
tant. Apart from its intrinsic interest, the investigation 
of persistence is further motivated by the observation, to 
be reported below, that P(to,t) is a sensitive probe for 
hidden temporal correlations in the growth process. This 
could shed some light on recent claims of nonuniversal 
scaling in discrete models associated with (|J) |J- 

Numerical results. Simulations were carried out in one 
dimension, using the strong-coupling discretization of the 
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KPZ equation proposed by Newman and Swift (NS) ||] 
as well as the restricted solid-on-solid (RSOS) growth 
model of Kim and Kosterlitz (pi. The NS model is 
a discrete-time model with continuous height variables 
h(x,t) evolving in parallel according to 

h(x, t + 1) = max[h(x, t), h(x — 1, t), h(x + 1, i)] + ij(x, t), 

(5) 

where the noise variables rj(x, t) are drawn from a Gaus- 
sian distribution with unit variance; the effect of using 
other distributions will be discussed at the end of the 
paper. 
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FIG. 1. Persistence probabilities for the NS model with 
Gaussian noise. The dashed lines have the slopes —1.21, and 
— 1.61, respectively. Inset: Steady state persistence probabili- 
ties obtained for a rough surface as initial condition (to = 10 4 
time steps, starting from a flat surface at time zero). The 
dashed line has the slope —2/3. 

We first checked relation (0) for growth from a rough 

(s) 

initial condition. Using the NS model we find 0\ J = 

6 ( l ] = 0.66 ± 0.03 (see the inset of fig. 0), which is in 
good agreement with (|J) and the exact value [3 = 1/3 of 
the one-dimensional KPZ equation . This confirms the 
conjecture |||| that the validity of (3]) is not restricted 
to Gaussian processes. 

To determine the nontrivial exponents 9±, the inter- 
face was started in the flat state, h — 0, and the per- 
sistence was measured from to — 1, i.e. one time step 
in the NS model and one deposition attempt per site in 
the RSOS model. The system size in the NS model was 
L = 16384 and averages were taken over 1000 indepen- 
dent runs. The measurements shown in fig. ^ yield the 
estimates 9 + = 1.21 ± 0.06, 0_ = 1.61 ± 0.08 which were 
obtained from fits in the regions 100 < t < 6000 and 
600 < t < 6000, respectively. The error bars are based 
on the fluctuations of the running exponent. The fact 



that 9+ < 9- is made plausible by noting that the non- 
linear term (A/2)(V/i) 2 in (Q) describes lateral growth, 
i.e. it causes, for A > 0, humps on the surface to extend 
sideways ]l|,0|. This tends to push negative segments of 
the surface (where h < h) upwards and leads to a faster 
decay of P_ . Since A > in the NS model |l2| , this ar- 
gument explains the inequality. It also suggests a faster 
decay of the height probability distribution in the direc- 
tion of positive heights when A < 0, in agreement with 
numerical [|l3j and exact jl4| results on related models. 

For the RSOS model the persistence probability was 
averaged over 124 runs, using a system size L = 5 x 10 5 
and 10 4 growth attempts per site. A fit to the data 
for times 500 < t < 5000 yields 9 + = 1.67 ± 0.07 and 
= 1.15 ± 0.08, with errors estimated from the ob- 
served fluctuations in the running exponent; the RSOS 
model has A < jnj , hence 9+ > 9- . Within error bars 
the exponents agree with those of the NS model. The 
numbers reported in the abstract represent the average 
of the values obtained for the two models. 
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= 0R/2)(1 - t'/t) 2/ ' A of the 
10°,10\ 10 2 ,10 3 ,10 4 . The 



FIG. 2. 
early-time behavior 1 — a(t,t') 
autocorrelation functions for t 

dashed line is the curve 0.85-(l-£'/t) 2/3 . Inset: The late-time 
behavior a(t,t') oc (t/t 1 )^ 1 of the autocorrelation function for 
t' = 10°, 10 1 , 10 2 , 10 3 , 10 4 . The dashed line is (t/t')~ l . 

Autocorrelation function. For a quantitative under- 
standing of the persistence exponents we need to char- 
acterize the shape of the KPZ autocorrelation function 
(|l|). It is useful to introduce two related functions: The 
temporal height difference correlation function 



C(t, t') = {{h{x, t) - h{x, t') - h{t) + h{t')) 2 

and the normalized autocorrelation function 
A{t,t') 



a(t,t') 



W(t)W(t') 



[t/t'fAit'/t). 



(6) 



(7) 



2 



The height difference function can be shown to have the 
form @|§ 

C(t,t') = c(t'/t)\t - t'\ 213 where the amplitude c(t'/t) 
tends to constants C2 and £2 for t'/t — > and t' ft — > 1, 
respectively (we take t > t'). Their ratio R — 02/02 is a 
universal number, whose value is numerically estimated 
[§U to be R = 1.8 ±0.1. Inserting this into (0) it follows 
that o(t,f) w 1 - (i?/2)(l - t'/t) 2 / 3 for t'/t -> 1 (fig. |). 

For t 3> t' the normalized autocorrelation function de- 
cays as a power law, a(t, t') ~ (t' /t) x , where the autocor- 
relation exponent A = 1 (fig- §)■ This is a special 
case of the general identity 



A = (3 + d/z 



(8) 



relating A to the surface dimensionality d, the roughness 
exponent and the dynamic exponent z; for the one- 
dimensional KPZ equation z = 3/2 [Q. Eq.(B) is valid 
for the linear growth equations considered in [2|, and we 
conjecture that it holds generally for rough growing inter- 
faces. To see this, note that a(t,t') measures the overlap 
between the height configurations at times t and t'. The 
overlap is the product of two factors: The lateral overlap 
between domains in the d-dimensional substrate space, 
which is of the order (£CO/£C0) with the dynamic cor- 
relation length |l[] ~ t x l z , and the horizontal overlap 
W{t')/W{t) ~ {t'/tf. 

Analytic bounds. Having established the behavior of 
the autocorrelation function for short and long times, the 
method of Ref. Q can be used to obtain an estimate for 
the persistence exponents #kpz of the Gaussian stochas- 
tic process associated with the KPZ autocorrelation func- 
tion <ikpz(£j O- We first make the process stationary by 
passing to logarithmic time . Setting T = In t we ob- 
tain ciKPz(t,t') = /kpz(T — T'), where the function /kpz 
has the limiting behaviors 



/kpz(T) 



1 - {R/2)\T\ 2 / 3 

„-T 



T 
T 





00, 



(9) 



and the persistence probability decays exponentially, as 
e -0KPzT _ rp ne f unc ^j on f Kpz w iH be compared to the 

logarithmic time autocorrelation function /g(T) of the 
linear growth equation with roughness exponent /? = 1/3, 
which is given by 



fp(T) = cosh(T/2) 



2/3 



|sinh(T/2)| 



■2U 



(10) 



The value /3 = 1/3 can be realized e.g. by the linearized 
version of (Q) with spatially correlated noise. For short 
and long times fp behaves as 1 -2- 2/3 |T| 2/3 and e _ ( 1_ #> T , 
respectively. Thus, at least in these limiting regimes, 
we have /kpz < A/3- If this is true for all T, a com- 
parison theorem due to Slepian [ fl6[ states that the per- 
sistence probabilities of the Gaussian processes associ- 
ated with /kpz and A/3 satisfy the same inequality, and 
hence 0kpz > #1/3 • In @ it was shown that 9y / 3 > 



(2/3)2 1 / 2 w 0.942809..., while the numerical solution of 
the linear growth equation yields 1 / 3 = 0.95 ± 0.05. 

More precise bounds can be derived by determining 
optimal scaling parameters & max and 6 m i n such that the 
inequalities /i/ 3 (&m ax 7 1 ) < /kpzCH < fi/ 3 (b m mT) hold 
for all T ||; Slepian's theorem then implies & m i n #i/3 



'KPZ 



< 



maxt'l/3' 



Comparing (^) and ( fio| ) and using the 



numerical estimate R « 1.8 we obtain the necessary con- 
ditions 6 max > 2(i?/2) 3 / 2 w 1.7 and 6 min < 3/2. Inspec- 
tion of the numerically determined function okpz yields 



the optimal parameters b„ 



1.35 and b n 



4, hence 



we conclude 1.28 < #kpz < 3.8. The fact that our nu- 
merical estimate for 9 + (which dominates the total persis- 
tence probability P + + P_ for long times) weakly violates 
the lower of these bounds appears to be a manifestation 
of the non-Gaussian character of the KPZ interface fluc- 
tuations. 
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FIG. 3. Persistence probabilities in the NS model with 
symmetric discrete noise. Inset: Period of the first (circles) 
and second (squares) oscillation as a function of q c — q. Using 
three-parameter fits we obtain q c — 0.542 in the first and 
q c = 0.537 in the second case. The dashed and the dotted 
line have the slopes —1.355 and —1.488, respectively. 

Non-Gaussian noise. Newman and Swift J9| recently 
reported a remarkable dependence of the roughness ex- 
ponent of high-dimensional KPZ growth on the probabil- 
ity distribution p{rf) of the noise variable i] in <^). Since 
the persistence probability (in contrast e.g. to the surface 
width or the autocorrelation function) probes temporal 
correlations of arbitrary order |l(J , it seemed natural to 
ask whether similar effects could be identified in low di- 
mensions using P±(to,t) as a dynamic observable. 

In fig. |H we show measurements of the persistence 
probabilities for the NS model with a symmetric, dis- 
crete noise distribution p{rj) — (l/2)[5(?j — 1) + 5(rj + 1)]. 
Instead of the expected power law decay (g), the P± 
show distinct oscillations, indicative of an intrinsic time 
scale t w 100; we emphasize that no oscillations were 
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found in the surface width or other more conventional 
quantities. The origin of the time scale r can be ex- 
plained by considering the general class of distributions 

Pq(v) = Q^iv ~ 1) + (1 — l)d(v + 1) where < q < 1. 
With this noise distribution the NS model (||) becomes 
equivalent to the polynuclear growth (PNG) model stud- 
ied by Kertesz and Wolf |l7[] , which is known to undergo 
a directed percolation transition at q c ~ 0.539; for q > q c 
the surface does not roughen. The transition is associated 
with a diverging correlation time r ~ \q — q c \~ Vt , where 
i>t ~ 1.73 is the temporal correlation length exponent of 
directed percolation. Measurements of the period of os- 
cillations in P± as a function of q show a similar power 
law divergence, confirming the identification of r as the 
directed percolation correlation time (inset of fig. ||). 

We have found similar, though less pronounced oscilla- 
tions for the continuous, bimodal noise distributions ||] 
Pa(v) ~ (1 — \ r l\) a with a > — 1. Further work is required 
to clarify the role of temporal correlations (induced by 
the parallel updating scheme (||)) in giving rise to the 
non-universal exponents observed by NS. The problem 
may indeed become more severe in higher dimensions due 
to the decrease of the directed percolation threshold q c ; 
in 2+1 dimensions q c ~ 0.272 [l7|j , hence for the symmet- 
ric discrete distribution (q = 1/2) the NS model resides 
in the smooth phase. 

Acknowledgements. J.K. acknowledges useful discus- 
sions with Alan Bray and financial support by DFG 
within SFB 237 H.K. acknowledges support by the Ger- 
man Academic Exchange Service within the Hochschul- 
sonderprogramm III. 



[13] J.M. Kim, M.A. Moore and A.J. Bray, Phys. Rev. A 44, 
2345 (1991). 

[14] B. Derrida and J.L. Lebowitz, Phys. Rev. Lett. 80, 209 
(1998). 

[15] S.N. Majumdar and C. Sire, Phys. Rev. Lett. 77, 1420 
(1996); S.N. Majumdar, C. Sire, A.J. Bray and S.J. Cor- 
nell, ibid. 77, 2867 (1996); B. Derrida, V. Hakim and R. 
Zeitak, ibid. 77, 2871 (1996). 

[16] D. Slepian, Bell Systems Tech. J. 41, 463 (1962). 

[17] J. Kertesz and D.E. Wolf, Phys. Rev. Lett. 62, 2571 
(1989); C. Lehner, N. Rajewsky, D.E. Wolf and J. 
Kertesz, Physica A 164, 81 (1990). 



[1] J. Krug and H. Spohn, in: "Solids Far From Equilib- 
rium", C. Godreche, ed., Cambridge University Press, 
Cambridge (1991), pp. 479; T.J. Halpin-Healy and Y.C. 
Zhang, Phys. Rep. 254, 215 (1995). 

[2] J. Krug, H. Kallabis, S.N. Majumdar, S.J. Cornell, A.J. 
Bray and C. Sire, Phys. Rev. E 56, 2702 (1997). 

[3] J. Krug, Markov Proc. Rel. Fields (in press). 

[4] M. Kardar, G. Parisi and Y.C. Zhang, Phys. Rev. Lett. 
56, 889 (1986). 

[5] J. Krug, Phys. Rev. A 44, R801 (1991). 

[6] J. Krug, P. Meakin and T. Halpin-Healy, Phys. Rev. A 
45, 638 (1992). 

[7] M. Krech, Phys. Rev. E 55, 668 (1997). 

[8] H. Kallabis, Ph.D. Dissertation, Universitat Duisburg 
(1997). 

[9] T. J. Newman and M.R. Swift, Phys. Rev. Lett. 79, 2261 
(1997). 

[10] J.M. Kim and J.M. Kosterlitz, Phys. Rev. Lett. 62, 2289 
(1989). 

[11] J. Krug and H. Spohn, Phys. Rev. Lett. 64, 2232 (1990). 
[12] H. Kallabis, unpublished (1998). 



4 



